Effects of Disorder and Momentum Relaxation on the Intertube Transport of 
Incommensurate Carbon Nanotube Ropes and Multiwall Nanotubes 
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We study theoretically the electrical transport between aligned carbon nanotubes in nanotube 
ropes, and between shells in multiwall carbon nanotubes. We focus on transport between two 
metallic nanotubes (or shells) of different chiralities with mismatched Fermi momenta and incom- 
mensurate periodicities. We perform numerical calculations of the transport properties of such sys- 
tems within a tight-binding formalism. For clean (disorder-free) nanotubes the intertube transport 
is strongly suppressed as a result of momentum conservation. For clean nanotubes, the intertube 
transport is typically dominated by the loss of momentum conservation at the contacts. We discuss 
in detail the effects of disorder, which also breaks momentum conservation, and calculate the effects 
of localised scatterers of various types. We show that physically relevant disorder potentials lead 
to very dramatic enhancements of the intertube conductance. We show that recent experimental 
measurements of the intershell transport in multiwall nanotubes are consistent with our theoretical 
results for a model of short-ranged correlated disorder. 

PACS numbers: 73.63.Fg, 73.22.-f, 72.10.-d 
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I. INTRODUCTION 

The electronic properties of carbon nanotubes are the 
source of much very interesting physics. Not only are car- 
bon nanotubes highly one-dimensional conductors, but 
the dependence of the nanotube bandstructure on the 
the nanotube's "chirality" leads to a range of different 
bandstructures: both metallic and semiconducting be- 
haviours arise and, even within each of these, the pos- 
sible bandstructures vary considerably^ The nature of 
the bandstructure has important effects on the electronic 
transport properties, not least the large difference in the 
effects of disorder on metallic nanotubes as compared to 
(doped) semiconducting nanotubes. ^'^''^ 

Theory shows that the bandstructure also plays key 
roles in the intertube transport between (the outer 
shells of) aligned nanotubes in a ropCfS, and in the 
transport between the shells of a multiwall nanotube 
(MWNT)&L2i2iia [Henceforth, we shall discuss these two 
cases together by referring to the transport as inievtuhe 
transport, with the two "tubes" to be understood as 
single-wall nanotubes (SWNTs) representing either the 
outer shells of two neighbouring nanotubes in a nanotube 
rope, or two neighbouring shells in a MWNT.] In ex- 
perimental systems of this type, the two tubes typically 
have differing chiralities and therefore differing band- 
structures. Even when the two tubes are both metallic, 
the typical case poses a complicated situation: the Fermi 
momenta do not coincide, and their translational peri- 
ods are incommensurate. Consequently, as a result of 
momentum conservation, one expects a strong suppres- 
sion of intertube tunnelling^ as compared to the simplest 
case of two identical tubes (with matched Fermi mo- 
menta and commensurate periodicity). Indeed, theory 
has led to the counter-intuitive prediction that increas- 
ing the length over which the two shells of a MWNT are 



coupled will reduce the intershell conductanceii 

A knowledge of the intertube conductance in nan- 
otube ropes and intershell conductance in MWNTs is 
of fundamental importance for understanding the trans- 
port mechanisms in these systems. However, direct ex- 
perimental studies of these quantities cannot easily be 
achieved. Very elegant recent experimental work has 
succeeded in probing the intertube/shell transport. Ex- 
perimental studies on intentionally disordered nanotube 
ropesii have been been used to estimate the intertube 
tunnelling conductance. Non-local resistance measure- 
ments between contacts on the outer shells of MWNTs 
have been used to determine the intershell conductivity^ 
(The development of telescopically extended MWNTsi^ 
opens the way to future studies of intershell conductance. 
However, at present the transport properties appear not 
to be dominated by the intershell conduct ance>i^) 

As emphasised in Ref. [T3 . the experimental measure- 
ments of intershell transport in MWNTs cannot be ac- 
counted for within existing theories. The experimen- 
tal systems have some degree of disorder, and existing 
theories of incommensurate tubes do not take disorder 
into account iSiL2i2ii£ Disorder is likely to have a very 
dramatic effect on the intertube/shell transport: disor- 
der breaks translational symmetry and therefore relaxes 
the momentum conservation which is responsible for the 
suppression^ of intertube transport. 

In this paper, we report calculations of the effects of 
disorder on the intertube conductance in carbon nan- 
otube ropes and MWNTs. In the rope geometry, we con- 
sider two nearest-neighbour tubes and study the trans- 
port between the outer shells of these tubes. In the 
multiwall geometry (in which the shells are coaxial) we 
study the transport between two neighbouring shells of 
the tube [thus we effectively consider a double wall nan- 
otube (DWNT)]. We focus on two metallic shells of arm- 



2 



chair and zigzag type, which provides an example of 
the typical case involving mismatched Fermi momenta 
and incommensurate periodicities. We show that de- 
fects lead to very dramatic enhancements of the inter- 
tube conductance. We quantify the effects of defects of 
different types, and discuss the results in connection with 
the experimental results of Refs.^ and ^3- Our results 
show that disorder has a large effect for both the inter- 
tube transport in disordered ropes measured in Ref. ITll 
and for the intershell transport in MWNTs measured in 
Ref.lll By comparing our theory to the experimentally 
determined intershell conductance of MWNTsiS we put 
constraints on the nature of the disorder in these sam- 
ples, and find parameter regimes of disorder for which our 
theoretical results are consistent with the experimental 
measurements ^ 



The paper is organised as follows. In Sj^we describe 
the model that we use, and discuss the importance of 
translational and rotational symmetries in the intertube 
coupling in aligned carbon nanotubes. We introduce a 
classification of the coupling strength between two shells 
of a MWNT based on the rotational symmetry. In ijllll we 
describe the Green's function approach used in our nu- 
merical studies. In ill VI we apply this to the case of clean 
(disorder- free) nanotubes. We demonstrate the impor- 
tance of momentum conservation for clean nanotubes. 
We show that momentum-relaxation at the contacts is 
the dominant effect controlling the intertube transport 
in clean nanotubes. The importance of our classifica- 
tion of the rotational symmetry in MWNTs is demon- 
strated through calculations of transport between nan- 
otube shells in a telescopic junction geometry. In Jvlwe 
turn to study the effects of disorder on the intertube con- 
ductance. We calculate the effects of localised defects of 
different types, and compare our results to the experi- 
mental measurements of Refs. and Finally in tiVIl 
we summarise the main results of the paper. 



II. NANOTUBE BANDSTRUCTURES, 
TRANSLATIONAL AND ROTATIONAL 
SYMMETRIES 

A. Theoretical Model 

We describe the bandstructure of the nanotubes using 
a one-band tight binding Hamiltonian, 

where Cia and c^^ are the annihilation and creation op- 
erators for an electron with spin polarisation a = ±1/2 
on the atomic site i. The atoms lie on a (rolled up) 
honeycomb lattice, with lattice constant a — 2.46A and 
hopping occurring only between nearest neighbour sites 
For quantitative estimates, we use t = 2.77eV con- 
sistent with the value fitted to the results of tunnelling 
experiments^ 

We consider intertube conductance between two car- 
bon nanotubes due to the addition of a tunnelling Hamil- 
tonian 

= - XI *y (cLcjT + c]aCm) (2) 

where i and j now label atomic sites on two different 
tubes. Following Ref. 0, we model the intertube hopping 
(between atoms at positions and rj) by 

t., =t±c-l'^--'-^l/^ (3) 

with tj_ = 492eV and S — O.SA. We have also performed 
calculations (not reported here) using a more sophisti- 
cated model^^ which accounts for the angular dependence 
of the hopping through a and tt bonds. We find similar 
results using either parameterisation in all cases we have 
studied (with realistic geometries). For the purposes of 
this paper, we therefore restrict attention to the results 
of the simpler model Q. 

In all the calculations we report, we treat the electrons 
as non-interacting. The results are thus valid at high 
temperatures, or for systems with many occupied modes, 
as in MWNTs, where the effects of interactions are sup- 
pressed. We discuss briefly the influence of interactions 
in El 



While preparing this work for publication, we learned 
of related worki4 where the importance of momentum re- 
laxation at the contacts of clean DWNTs was also recog- 
nised. Since our discussion of this aspect in ijlVI is some- 
what different from Ref. and since an understanding 
of the clean case is required in order to introduce our 
calculations of the effects of defects, we retain this dis- 
cussion. Our results in this section are consistent with 
those of Ref. B 



B. Bandstructures of Armchair and Zigzag 
Nanotubes 

It is well-known that the electronic properties of car- 
bon nanotubes depend on the chiral vector (ni,n2), 
which sets the quantised subbands and thus the 
bandstructurCfi In this work we focus on metallic arm- 
chair and zigzag nanotubes. 

Armchair tubes are described by {na,na) chiral vec- 
tors. The wavevector around the tube circumference 
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FIG. 1: (Color online) Schematic diagrams of the lattice 
structures (the tube axes run horizontally) and bandstruc- 
tures of the metallic bands for armchair and zigzag nanotubes. 
The repeat distances along the tube axes differ by a factor of 
\/3, so the tubes are incommensurate. Intertube tunnelling 
will be suppressed for Fermi energies close to E = 0, owing 
to the lack of energy-momentum conserving processes. 



takes the values k± = {2n/^/3a){q/na) with q = 
l,2,...2nQ. For each q there are two bands, labelled 
by ±, which have the energy dispersions 



1 + 4 cos 



(77r\ / ka\ . .-, f ka 
— cos — -)-4cos^ — 

(4) 

where k is the wavevector along the tube, with — tt < 
ka < TT. The ± bands with q = n are metallic, crossing 
E = Oa.tk = ±27r/(3a). 

Zigzag tubes are described by (^2,0) chiral vectors. 
The wavevector around the tube circumference takes val- 
ues k± = (2TT/a)(q/nz), with q = I, ...,2nz, and the sub- 
bands have dispersion relations 



C(fc) = ±t 



( ^/ika\ ( q-K 
4 cos — - — cos — 
2 i \n. 



4 cos 



q-K 
(5) 

with -tt/VS <ka< 7r/\/3. If Uz is an integer multiple 
of 3, the zigzag tube is metallic. The metallic subbands 
have q — 2n^/3, 4^2/3 (the ± bands correspond to right- 
and left-moving modes of these two subbands). The two 
metallic subbands both cross = at fc = 0. 

The lattice structures and metallic bands of armchair 
and zigzag nanotubes are illustrated in Fig. ^ 



C. Intertube Conductance: Momentum 
Conservation and Incommensurability 

We explore the effects of momentum (wavevector) 
conservation in the intertube conductance between two 
aligned nanotubes. First, in order to explain the impor- 
tance of momentum conservation, we consider the tubes 



to be sufficiently weakly coupled that this coupling can 
be described by perturbation theory. This is valid pro- 
vided the intertube conductance is small compared to 
/h. (Our numerical results will not be restricted to 
this perturbative limit.) Within perturbation theory, for 
transport between a Fermi point fcA on one tube and a 
Fermi point fee on the other there is a contribution to 
the intertube conductance of 



GkA,kB 



2e2 



r2 



H^vavb 27r 



\tkA 



where 



(feBlHrlfcA) 



(6) 



(7) 



is the tunnelling matrix element, va and wb are the ve- 
locities at the two Fermi points, and Ljv is the length 
over which the states |A;) are normalised. The total con- 
ductance is obtained by the sum of © over all initial and 
final Fermi points (fcA,fcB) associated with the different 
occupied subbands [we have suppressed band indices in 
Eqns.ElandCI. 

Momentum conservation arises from the symmetry of 
the tunnelling Hamiltonian under translation along 
the tube.^^ If the Fermi momenta of the two tubes are 
the same, (fcA = fce), then electrons at the Fermi energy 
can tunnel between the tubes conserving both energy and 
momentum: one expects a large intertube conductance. 
(For two tubes of the same type, the Fermi momenta 
coincide at all energies so intertube coupling is always 
strong.) If the Fermi momenta differ (kp^ 7^ fce), tun- 
nelling will be suppressed. In general, tubes of different 
chiralities have different dispersion relations and there- 
fore differing Fermi momenta. For example, as shown in 
Fig. ^ the Fermi momenta of the armchair and zigzag 
tubes differ in a range of Fermi energies around E — Q. 
(At Ep = the Fermi momenta are fc|f = ±27r/(3a) and 
fc|? = for armchair and zigzag tubes respectively.) 

Although one might conclude that if the Fermi mo- 
menta differ, intertube tunnelling will vanish identically, 
this is not always true. An important distinction arises 
between the cases in which the tubes have commensu- 
rate or incommensurate lattice structures. Denoting the 
repeat distances in the two tubes by a a and as, we say 
that the carbon nanotubes are commensurate if ap^/a^ 
is rational fraction, and that they are incommensurate if 
this is an irrational fraction. 

If the carbon nanotubes are commensurate, then the 
system has translational symmetry along its length. The 
minimum distance along the tubes under which the sys- 
tem is translationally invariant sets the size of the first 
Brillouin zone. If the Fermi momenta of the two tubes do 
not coincide within this first Brillouin zone, then there 
can be no tunnelling between the tubes that conserves en- 
ergy and (crystal) momentum. The intertube tunnelling 
will vanish. 

If the carbon nanotubes are incommensurate, then 
there is no finite distance under which the system is 
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translationally invariant: the size of the first Brillouin 
zone is vanishingly small, and the eigenstates cannot be 
characterised by a conserved crystal momentum. Never- 
theless, for weak coupling between the tubes, the Bloch 
states of the individual tubes still provide a useful start- 
ing point. Imagine that, in the first Brillouin zones of the 
individual tubes, the Fermi momenta of the two tubes do 
not coincide. Because the reciprocal lattice vectors are 
incommensurate, by adding sufficient reciprocal lattice 
vectors to each of the tubes one will be able to bring 
the Fermi momenta arbitrarily close. (That is, in an ex- 



tended zone scheme the Fermi momenta can become arbi- 
trarily close.) However, as shown by Maarouf et al^, the 
larger the number of reciprocal lattice vectors that must 
be added to achieve momentum matching, the weaker is 
the tunnelling between the two tubes. Ref. 3 provided an 
approximate calculation of the matrix element of the in- 
tertube coupling 121 between momentum eigenstates on 
two aligned nanotubes placed side-by-side (as in a carbon 
nanotube rope). Following this approach;^ which is valid 
in the limit that tunnelling is dominated by the points of 
closest approach, we find that the matrix element is 



t 



/ 27r 6^5 



(8) 



where Acc\i is the area of a unit cell, fcj^' are the 
(subband) wavevectors around the tube circumference, 
= + aC^ and k^ = k^ + [3G^ are the wavevec- 
tors along the tubes shifted by integer numbers (a, /?) 
of reciprocal lattice vectors along the tube (G'^,G^), 

= I I + I ^? I ' ^^"^ ^ distance between tubes 

at the point of closest approach. For larger numbers of 
reciprocal lattice vectors, increases and the tunnelling 
matrix element is exponentially suppressed. One expects 
the tunnelling to be dominated by processes involving the 
smallest number of additional reciprocal lattice vectorsi^ 
Two tubes of different chiralities typically have Fermi 
momenta that do not coincide, and typically are incom- 
mensurate. The armchair and zigzag tubes are an ex- 
ample of this typical case. Here, the smallest reciprocal 
lattice vectors are = 27r/a, G^'" = 27r/V3a. The 
condition for energy-momentum conservation is 

e^l [kf + aG^^) = e'l {k'f + /3G^") = Ep (9) 

where k^p^^^ are the Fermi momenta along the respective 
tube axes, and a and /3 are integers. Focusing on the 
dominant tunnelling process involving just one reciprocal 
lattice vector of each tube, we look for a solutions with 
a = ±1 and /? = ±1. Using the dispersion relations 
(|4I5|) we find ka — ±3.914 at energies Ep — zLEq where 
Eq = 0.247i. We shall show the importance of this energy 
in the numerical results in illVBI 



D. Multiwall Geometry: Rotational Symmetry 

In the case of MWNTs, intertube conductance is 
strongly affected not only by momentum conservation 
along the tube, but also by angular momentum conserva- 
tion around the axis of the tube. This leads to additional 
selection rules in the tunnelling matrix elements which 
depend sensitively on the chiralities of the two tubesn^ii^ 



We discuss the implications of this symmetry for the 
tunnelling matrix element |(7J) between armchair and 
zigzag shells. The subband wavefunctions vary as e'*^^^, 
where y is the coordinate around the circumference 
of the tube. For an {na,na) armchair tube, k± — 
271 /(^/3a){qa/na), and the lattice is invariant under a 
translation by \/3a in y. Using 9f — yi/ {s/iana) = i/ria 
to label the angular position of the unit cell i (i — 
1, . . .Ua), then the subband wavefunction is simply e"^"^» . 
Similarly for the zigzag tube, 0| = yj /{ariz) — 3 jn^^ and 
the subband wavefunction is e''-^j for the unit cell la- 
belled by J = 1, . . .n2. Therefore, the tunnelling matrix 
element contains an angular contribution of the form 

(*.|Ht l^a) (X ^t(0,'' - ^|)e'«"«"e-'«^«.' 

= ^ ^t„e-"(^?-«.^)e'«"«?e-'«^«.^ 

where t(6'f — ^|) is the function obtained by expressing 
Q in terms of the angular separation between atoms in 
the two shells, which in the last line has been expressed 
as a Fourier series due to its periodicity. Performing the 
sums over i and j, one finds that the contribution of the 
harmonic t„ is nonzero only if 

Condition A: mod(n — (/a, na) = (10) 
Condition B: mod(n — (72, n^) = (11) 

If both of these conditions are met, then modes qa and 
are connected by a non-zero tunnelling matrix element. 
Restricting attention, now, to the metallic bands in both 
tubes - that is, qa — na and q^ = 2nz/3 or q^ = 47X2/8 - 
the conditions (|10|l . (|11() may be rewritten 



A: n = aua 
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FIG. 2: Geometry of the system that we study, illustrating the 
labelling of the leads and the length L of the central section 
where intertube tunnelling can occur. The length Lm, which 
includes the leads, is used for normalisation and is assumed 
large. 



We classify the nature of the rotational symmetry of 
the tubes in three cases. Which case applies depends on 
the ratio na/riz — p/q where p and q are natural numbers 
with no common factor. 

(i) "Zero": No value of n that satisfies condition B 
satisfies A. This is the case if q is not divisible by 3. The 
tunnelling matrix elements between the metallic bands 
of the armchair and zigzag tubes vanish identically. 

(ii) "Strong": All values of n that satisfy condition 
B, satisfy A. This is the case if q is divisible by 3 and if 
p = 1. The metallic bands of the armchair and zigzag 
tubes couple strongly, containing contributions from all 
possible angular harmonics of the tunnelling operator. 

(iii) "Intermediate": Some values of n that satisfy 
condition B satisfy A. This is the case if q is divisible by 
3 and if p 7^ 1. The metallic bands of the armchair and 
zigzag tubes are coupled, but contain contributions from 
only a subset of the angular harmonics of the tunnelling 
operator. 

That rotational symmetry can lead to vanishing matrix 
elements (the "zero" case) is well-knownii^ However, as 
we show below in our numerical results, there can be 
very substantial quantitative differences between cases 
we discriminate as "intermediate" or "strong" coupling. 



III. CONDUCTANCE USING GREEN'S 
FUNCTIONS 

In our numerical studies, we calculate the conductance 
of the double nanotube system within the Landauer- 
Biittiker framework, with the scattering matrix calcu- 
lated numerically using a Green's function approachiiS 
The geometry of the system studied is shown in Fig. |21 
It consists of two parallel nanotubes, each of which we 
shall consider to be either armchair or zigzag type. The 
geometry allows a description of (the outer shells of) two 



neighbouring tubes in a rope geometry (if the tube axes 
do not coincide) , or two neighbouring shells in a MWNT 
(if the axes do coincide). 

We calculate the conductance using a scattering ap- 
proach. For this purpose, we consider leads to be at- 
tached to the ends of each nanotube. The leads are ideal, 
consisting of a semi-infinite section of clean nanotube, 
and there is no intertube coupling between the leads. 
The intertube coupling is confined to the central section 
of length L. In this region the system has Hamiltonian 
Tic, which describes atomic sites that experience intra- 
and inter-tube hopping and where there can be diago- 
nal or off diagonal disorder. The nature of the intertube 
scattering in this region determines if we are considering 
the rope or multiwall geometry. 

Following the approach described in Ref. 19, the over- 
all Green's function at energy E can be obtained via a 
matrix inversion 



Gc = [EI -Hc-n 



(12) 



where the matrix has dimension C x C, with C the num- 
ber of sites in the central region. [The typical system 
sizes we study have on the order of C = 9000 atoms. 
Matrix inversion is done using standard LaPack code.] 
The effects of the leads are taken into account via a self- 
energy term S = YipT^g^Tp which can be viewed as an ef- 
fective Hamiltonian describing the coupling between the 
conductor and the four leads (labelled by p = 1, ... 4). In 
this expression, gp is the end Green's function of lead p 
and Tp is a matrix describing the coupling of this lead to 
the central region. Although the leads are semi-infinite, 
the end Green's function has finite dimension, set by the 
number of sites on the end of the lead. 

We restrict attention to the case in which the energy 
is sufficiently close to = that there are only two con- 
ducting modes in each of the metallic nanotube leads. 
This allows a description of transport properties when 
the energy E is less than the gap to the next subband. 
In this case, it is sufficient to calculate the Green's func- 
tion for the leads including only these metallic modes. 
We have determined the end Green's function analyti- 
cally for the metallic modes of the armchair and zigzag 
tubes (the calculation is described in Appendix IVIip . al- 
lowing an analytic determination of the self-energy enter- 
ing Eqn.^l Note that in the "central" region the system 
is solved fully (numerically), including all subbands. 

For the range of energy that we study, the transmis- 
sion matrix is an 8 x 8 matrix, T^ p.^^i pi, describing the 
transmission from two incoming modes (labelled by /i) 
in each of the four leads (p = 1, 2, 3, 4) to the 2x4 out- 
going modes (with mode ^' and lead p') The nature of 
the mode depends on whether the lead is of armchair or 
zigzag type: for an (jiaTna) armchair tube, /i labels the 
± bands with qa = ria] for a (0, n^) zigzag tube, fi labels 
qz = 2n^/3,4nz/3. We obtain the scattering matrix for 
currents, following the approach described in Ref. for 
a continuum system applied to this discrete case. As in 
the case of the continuum, it is important to convert the 
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matrix for scattering amplitudes into the matrix for scat- 
tering fluxes using the velocities of the initial and final 
modes. From the scattering amplitude for the flux, one 
obtains the transmission probabilities, the final results 
for which are given in ADpendi? iVIIII 

Finally, using the Landauer-Biittiker approach, the 
conductance is obtained from the sum of the transmis- 
sion probabilities. For example, the intcrtube conduc- 
tance between tube-A {p = 1,2) and tube-B {p = 3,4) 
as measured in an experiment with contacts attached to 
all four leads, involves the total transmission coefficient 
between all modes in the tube-A and all modes in tube-B 

/i;p=l,2 ii' ;p'=3,4 

The prefactor of 2 arises from spin degeneracy. The max- 
imum intertube conductance in this case is 8e^/h. 
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FIG. 3: (Color online) Two-terminal conductance of a (5, 5) 
armchair nanotube with a single vacancy (solid line) and a 
(9,0) nanotube with a single vacancy (dashed line). The tube 
length [L = 31a) used in these calculations is sufficiently long 
that the results are characteristic of the bulk behaviour. 



A. Application: Conductance of a SWNT with a 
Vacancy 

As a first application of our numerical method, we 
study the two-terminal conductance of single SWNTs 
(both armchair and zigzag) in the presence of a vacancy. 
This is achieved simply by turning off intertube tun- 
nelling, t± = 0, and calculating the scattering between 
the two ends of the individual tubes. Thus, the two- 
terminal conductance of tube-A is 

(Jintra — / ^ J/j,l;;i',2 

which has maximal value Ae^ /h for a clean metallic nan- 
otube (two propagating modes and twofold spin degen- 
eracy) . 

Related calculations of the effects of vacancies on arm- 
chair SWNTs have been reported previouslyi22*2iiS242^ As 
in those studies we do not include relaxation of the lat- 
tice around the vacancy, but leave the honeycomb lattice 
undistorted. Following Ref. [23 we model the vacancy by 
imposing a huge potential on the vacancy site (lO^eV), 
which was shownS^ to be equivalent to explicit discon- 
nection of the bonds to this site. 

In Fig. 131 we present results for the conductance of sin- 
gle metallic armchair and zigzag carbon nanotubes in the 
presence of a single vacancy. Our calculations reproduce 
the results previously found by Igami et al?^ for arm- 
chair tubes. In addition, we present results for a metal- 
lic zigzag nanotube. The results in both cases (armchair 
and zigzag) are very similar: the conductance drops from 
close to its ideal value of 4e^//i at the edge of the second 
conduction band, to 2e^//i at E = 0.— 

One can understand that the minimum value of the 
conductance at = is Gi„tra = 2e^//i by noting that, 
in both armchair and zigzag cases, a linear combination of 



the two degenerate right-moving (or left-moving) states 
of an ideal (defect free) nanotube can be formed that has 
zero amplitude on a subset of the lattice sites. Choosing 
the linear combination such that the amplitude is zero 
on the site containing the vacancy, this mode of the ideal 
tube is unaffected by the presence of a vacancy and there- 
fore contributes 2e^//i to the two-terminal conductance. 
That our numerical results show that the conductance 
falls to precisely this minimum value at E ~ Q indicates 
that the other (orthogonal) mode is completely reflected 
by the vacancy. For an armchair tube, the wavefunction 
of the transmitted mode is periodic along the length of 
the tube every three lattice constants [it is a superpo- 
sition of waves with k — -|-27r/(3a) and k = — 27r/(3a)]. 
Therefore, one expects that it should be possible to intro- 
duce several vacancies, at positions chosen judiciously to 
lie at the set of zeroes of this transmitted mode, without 
having any effect on its transmission. We have checked 
this expectation numerically by introducing several va- 
cancies in the carbon nanotube on sites respecting the pe- 
riodicity of the wavefunction: a,t E = Q the conductance 
remains 2e^//i, indicating that the transmitted wave in- 
deed does not feel the potentials on these sites. 

IV. INTERTUBE CONDUCTANCE BETWEEN 
TWO CLEAN SWNTS 

We now turn to consider the intertube conductance 
between aligned nanotubes in both rope and multiwall 
geometries. From the discussions of iJII CI we expect that 
momentum conservation and, in the case of MWNTs, 
the rotational symmetry will play important roles in de- 
termining the magnitude of the intertube conductance. 
Here, we investigate these effects, focusing on the case 
of clean tubes in the absence of defects. A discussion 
of the effects of defects will be given in We begin 
our discussion by considering the effects of the loss of 



7 



momentum conservation at the contacts. 



A. Contact Effects 

The presence of current contacts at the two ends of a 
carbon nanotube sets a finite total length of the system, 
L. In Ref. 7 it was recognised that as a consequence of 
this finite length, L, momentum conservation will be re- 
laxed on a wavevector scale of order 1/L. This led to the 
prediction^ of a decrease in the intershell conductance of 
a MWNT with increasing length L, as momentum con- 
servation becomes more accurately satisfied. 

Another effect is that in the vicinity of the contacts 
themselves the translational symmetry is broken. We 
denote the (minimum) intrinsic lengthscale of the con- 
tact Lc- this could be the size of the contact, or the 



lengthscale over which the electronic nanotube density 
changes between a region under the contact and a re- 
gion away from the contact. The variation of the nan- 
otube properties on will promote a loss of momen- 
tum conservation up to a wavevector scale of about 1 / L^. 
Since 1/Lc ^ this change in momentum will typ- 
ically be larger than that arising from the finite length 
of the tube. We find that momentum relaxation at the 
contacts typically dominates the transport properties be- 
tween clean nanotubes with mismatched Fermi momenta. 
The contact-induced contribution to the intcrtubc con- 
ductance was not considered in Ref. 0- 

Since we are interested in the general features of mo- 
mentum (non-) conservation at the contacts, we intro- 
duce a simple model that incorporates a contact length- 
scale Lc as well as the total length L. Specifically, we 
replace (O by 



J 



(14) 



where /(x) = l/(e^-|-l) and the distances along 

the two tubes. This describes intertube tunnelling over 
a total length L, but with the tunnelling rising smoothly 
over a lengthscale Lc from zero (far inside the "leads" ) to 
its maximal value. Although this model for smooth con- 
tacts (|14|) does not accurately mimic experimental situa- 
tions, it provides a convenient way in which to control the 
extent of scattering at the contacts. We shall make use 
of this model (fnji throughout the paper, choosing suf- 
ficiently large Lc to eliminate the effects of momentum 
relaxation at the contacts and allow the bulk properties 
to be studied. 

The qualitative effects of L and Lc can be under- 
stood in simple terms by considering the perturbative 
expression for conductance © for a one-dimensional 
continuum model with tunnelling amplitude t{xi^X2) — 
t(x)5(xi — X2), with 



t{x) = io [e (x) - e (x - L)] * {x) 
where the star denotes convolution, and 

1 2e"^/-^'^ 



(1 -f- e-^/L^y 



(15) 



(16) 



For L ^ Lc and S <ti Lc this is an accurate approxima- 
tion to the functional form (|14|l . The matrix element of 
the tunnel coupling from states /ca to fee is 

ifcAfcB = Y^J I t{xi,x2)e-"'^^'e"'^^'dx,dx2{17) 



Ak 



(18) 



where Ak 



kx and 

7rAfcLc(l + iAkLc) 



sinh {irAkLc 



(19) 



is the Fourier transform of gL^{x), which was obtained 
by closing the contour in the correct half-plane and sum- 
ming the residues of an infinite number of third order 
poles at a; = {2n + l)i7r. This matrix element enters the 
perturbative expression for conductance © to give 



HHavb {Akf 
In the case of momentum conservation, Ak 



(20) 



0, 



G 



'-0 



2^3 



vavb 



(21) 



That is, the conductance grows like the square of the 
tunnelling length L, and is independent of the contact 
length Lc- these results are consistent with the conduc- 
tance arising from the bulk of the system, unaffected by 
the contact regions. 

For the situation of momentum non-conservation, 
Afc 7^ the conductance (|20|) is an oscillatory function 
of the total length of the system, L. This contribution 
arises from tunnelling at the contacts, with the oscilla- 
tory dependence a consequence of interference between 
the tunnelling at the contacts at two ends of the sam- 
ple. Consistent with this explanation, the amplitude of 
this oscillatory contribution to conductance depends sen- 
sitively on the contact lengthscale Lc, with exponential 
suppression for large AkLc- 
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FIG. 4: Intertube conductance between (10,10) and (18,0) 
tubes, in the rope geometry with an intertube distance of 6 = 
3.4A. The total length of the carbon nanotubes is L = 87a, 
and intertube tunnelling is introduced smoothly from zero 
over a length Lc = 2a. The peaks at E/t = ±0.25 correspond 
to energy-momentum conserving processes described in the 
text. 



B. Rope Geometry 

We start by discussing the intertube conductance be- 
tween the outer shells of two neighbouring nanotubes 
that lie side by side in a nanotube rope. We choose an 
intertube separation oi b — 3.4A which is typical for a 
nanotube rope. 

If the two tubes are identical, they are commensurate 
and momentum conservation is satisfied at all energies. 
We therefore expect strong coupling. Indeed, we find in 
our numerical calculations (not shown) that for two arm- 
chair tubes or two zigzag tubes the intertube conductance 
is large, reaching of the order of e^//i even for very short 
tunnelling regions between the tubes. 

Typically, two tubes of different chirality have mis- 
matched Fermi momenta and are incommensurate. As 
discussed in ijll CI we expect tunnelling to be strongly 
suppressed. We investigate this effect for the case of 
(10,10) armchair and (18,0) zigzag tubes. In Fig^ 
we show the intertube conductance as a function of 
the Fermi energy E. (We have suppressed the effects 
of momentum relaxation at the contacts by choosing 
Lc — 2a.) The conductance is very small for most ener- 
gies. However, there are two peaks (with sidebands), at 
E « ±0.25t. At these energies tunnelling can occur by 
the momentum-conserving tunnelling process discussed 
at the end of ^11 Cl in which a single reciprocal lattice vec- 
tor is added to the Fermi momenta of both armchair and 
zigzag tubes. The presence and positions of these peaks 
are consistent with the theory of Rcf. 5 and formula © 
which showed that the dominant tunnelling processes in- 
volve the smallest number of reciprocal lattice vectors. 
However, the heights are not well reproduced by ((Sj). In 
particular, the large asymmetry in the conductance at 
the peaks E = ^Eq, by several orders of magnitude, is 
not reproduced by Eqn. |S1 These failings of Eqn.|Slshow 
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FIG. 5: (Color online) Length dependence of the intertube 
conductance for (10, 10) and (18, 0) nanotube in a rope geom- 
etry a distance b = 3.4A apart. Results are shown as a func- 
tion of the total length, at E = —Eo,0,+Eo (solid, dotted, 
dashed, respectively). The main figure, which plots G/{La)^ , 
indicates clearly that G ^ for E = Eq. The inset, which 
plots G, indicates that G ^ also for E = —Eg but not for 
_B = 0. These results are consistent with a momentum con- 
serving tunnelling processes a,t E — ±Eo but none at i5 = 0. 
(A contact length of Lc = 2a has been chosen to minimise 
effects of momentum relaxation at the contacts.) 



that one must go beyond this approximate expression to 
achieve quantitatively accurate values for the intertube 
conductance. 



The continuum model of i|IV Al accounts for the qual- 
itative features of the length-dependence of the conduc- 
tance of the (10, 10) to (18,0) system. As illustrated in 
Fig. \^ a.t E = ±0.25t we find that G cx L'^, charac- 
teristic of a bulk momentum-conserving tunnelling pro- 
cess [see Eqn. At E — the conductance oscillates 
with length, indicating that no momentum-conserving 
processes are contributing. The observed period of the 
oscillation is consistent with the form sin^(AfcL/2), with 
the change in wavevector \Ak\ = (2 — •\/3)27r/a expected 
from the difference between the Fermi wavevectors of the 
armchair and zigzag tubes up to a shift by reciprocal 
lattice vectors. Furthermore, the contact-induced con- 
tribution to the conductance is a strong function of La- 
in Fig. we show this dependence, plotting the conduc- 
tance at i? = (averaged over two periods of oscillation) 
as a function of contact length. As shown in Fig. the 
rapid decay of conductance with Lc can be fit to the 
form given by Eqn. 1201 Deviations from this form indi- 
cate the breakdown of the continuum approximation for 
small Lc S, «• 
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FIG. 6: (Color online) Mean intertube conductance at iJ = 
between (10,10) and (18,0) carbon nanotubes in the rope 
geometry a distance b = 3.4A apart, as a function of the 
contact length Lc- The strong dependence on Lc indicates 
that the conductance is dominated by momentum relaxation 
at the contacts. The fit (dashed line) is to the functional form 
of Eqn. I2S1 




FIG. 7: (Color online) Mean intertube conductance at _E = 
between (7, 7) and (21, 0) carbon nanotubes in a multiwall ge- 
ometry as a function of the contact length Lc- The strong de- 
pendence on Lc indicates that the conductance is dominated 
by momentum relaxation at the contacts. The fit (dashed 
line) is to the functional form of Eqn. 1201 



C. Multiwall Geometry 

For the multiwall geometry, in addition to momentum 
conservation along the tube, the classification of the ro- 
tational symmetry, 311 Dl plays a key role. In order to 
illustrate this we have calculated the intertube conduc- 
tance (|13|l for clean DWNTs of the three different kinds: 
"zero" (3,3) in (15,0); "intermediate" (5,5) in (18,0); and 
"strong" (7,7) in (21,0). These are chosen to have similar 
intershell spacing (3.8A, 3.7A and 3.5A respectively). In 
each case, since the tubes have mismatched Fermi mo- 
menta, the conductance is small owing to momentum 
conservation along the tube. However, the relative sizes 
of the conductances in these three cases differ greatly: 
for the "zero" case, we find vanishing intertube con- 
ductance to machine accuracy; the conductance of the 
"strong" coupling case is larger than that of the "inter- 
mediate" case by about 8 orders of magnitude (results 
are not shown here). This very significant difference in 
conductance indicates the importance of distinction be- 
tween "strong" and "intermediate" coupling cases. 

As in the rope geometry, the conductance is a strong 
function of the contact length Lc- In Fig|7|this depen- 
dence is shown for the average conductance at i? = of 
the (7,7) to (21,0) DWNT (which is the "strong" cou- 
pling case in terms of the rotational symmetry). Note 
that, for small values of the contact lengthscale, a con- 
ductance of order 0.05e^/h can be induced. This is a 
sizable conductance, given that it is due only to currents 
that flow in the vicinity of the contacts. Our results 
emphasise that to understand the intertube conductance 
between clean tubes, with mismatched Fermi momenta 
and incommensurate periodicity, requires a very careful 
study of the contacts: the conductance is dominated by 
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FIG. 8: Schematic picture of a telescopic junction. 

the current flowing in the vicinity of the contacts. 



D. Telescopic Junction 

A particularly simple contact geometry between two 
neighbouring shells of a MWNT or DWNT can be formed 
by telescopic extension of the inner tubes^ and making 
separate contacts to the inner and outer shellsi2& Current 
that is passed between these two contacts is then forced 
to flow between the two shells. (See Fig|Sl) This geom- 
etry involves a very abrupt change in the properties of 
the tubes at the free ends in the inner and outer tubes. 
We therefore anticipate that the loss of momentum con- 
servation at these edges will play a significant role in the 
intertube conductance. 

We have evaluated the conductance for telescopic junc- 
tions formed from clean DWNTs. To model this situa- 
tion within our approach, the leads p = 2, 3 in Fig[21 
are disconnected (the r matrices discussed in illlll are 
set to zero) and the conductance between leads 1 and 4 
is calculated. Similar studies have been performed pre- 
viously for certain commensurate and incommensurate 
tubes, 
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FIG. 9: Conductance of a telescopic junction between two 
incommensurate (6, 6) to (18, 0) carbon nanotubes, as a func- 
tion of the overlap length, £o. By our classification of the 
rotational symmetry the coupling is "strong" . The two tubes 
are strongly coupled even for small overlap lengths. 



Here we focus on the case of incommensurate tubes 
with mismatched Fermi momenta, and discuss the influ- 
ence of the rotational symmetry classification of ijllDI A 
telescopic junction does not break the rotational symme- 
try, so this classification remains relevant. We present 
results for both "intermediate" [(10,10) to (9,0)] and 
"strong" coupling [(6,6) to (18,0)] in Figs. and UHl 
(For the "zero" case, we find that the conductance van- 
ishes, as expected by the symmetry analysis.) For the 
"strong" coupling. Fig. we indeed find very efficient 
coupling. For even small overlap lengths of a few lattice 
constants, the conductance is of order e^/h. Attempting 
to describe the system as two standing waves in different 
tubes is no longer valid in this regime: since the inter- 
shell conductance is of order e^/h the tubes are strongly 
coupled and perturbation theory is inaccurate. For inter- 
mediate coupling, Fig llOl we find an oscillatory intershell 
conductance but of a very much smaller magnitude. In- 
creasing the overlap does not increase the conductance, 
consistent with the fact that this is a system with mis- 
matched Fermi momenta, for which there are no bulk 
momentum-conserving tunnelling processes. These re- 
sults again highlight the striking difference between ro- 
tational symmetry classification as "strong" and "inter- 
mediate" . In the former case, the relaxation of momen- 
tum at the junction between SWNT and DWNT region 
is enough to cause a conductance of order e^/ft; in the 
latter, the intershell conductance is many orders of mag- 
nitude smaller. 



V. EFFECTS OF DEFECTS ON INTERTUBE 
CONDUCTANCE 

In ijlVl the intertube conductance was studied for clean 
(defect-free) tubes. Momentum conservation was shown 
to play a very important role in the intertube conduc- 



FIG. 10; Conductance of incommensurate (10,10) to (9,0) 
carbon nanotubes for a telescopic junction, as a function of 
the overlap length £o. By our classification of the rotational 
symmetry the coupling is "intermediate" . The two tubes re- 
main weakly coupled even for large overlap lengths. 

tance between tubes with mismatched Fermi momenta. 
In real samples, electrical conductivity typically depends 
vitally on the number and type of defects present. One 
can expect defects to act as scattering centres and to 
relax momentum conservation. This general point was 
noted in Ref. ,5,. However, we are not aware of any stud- 
ies taking into account the effects of defects on intertube 
or intershell transport. In this section we provide cal- 
culations of the effects of defects of different types on 
the conductance of incommensurate nanotube ropes and 
MWNTs. In particular, we discuss the effects of individ- 
ual scatterers (vacancies and other localised defects) on 
intertube transport in cases where the clean tubes would 
have small intertube conductance. We shall then discuss 
how these results can be extended to a finite density of 
scatterers. 

In all the studies we report below, we use a contact 
length of Lc = 2a to suppress the effect of momentum 
relaxation at the contacts and the total length of the 
central region is made sufficiently long (of order L = 
100a) that we are able to extract the inffuence of the 
impurity in the bulk system without finite size effects. 



A. Vacancies 

While rare in natively grown nanotubes, vacancies can 
be introduced by sputteringiii Here we determine the 
effects of vacancies on intertube conductance in ropes 
and MWNTs. 



1. Rope Geometry 

In Fig^Jwe present results for the effects of a vacancy 
on the conductance between (10,10) and (18,0) carbon 
nanotubes in the rope geometry. The nanotubes both 
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FIG. 11: (Color online) Intertube conductance between a 
(10,10) armchair and (18,0) zigzag in the rope geometry a 
distance b = 3.4A apart. Results are shown both for clean 
carbon nanotubes (dashed line) and with a single vacancy 
placed on a site halfway along the armchair tube (solid line). 
The inset shows the dependence of the conductance at _E = 
on the location of the vacancy at the atomic sites around the 
circumference of the armchair tube. 

have diameters of about 1.4nm and are placed 3.4A apart. 
The vacancy, which we introduce in the same way as in 
^III Al is placed on the armchair tube at the point of 
closest approach. 

In Fig. ^2 the dashed line shows the intertube con- 
ductance for clean tubes, which was discussed in ijlV Bl 
and shown in Fig^ (on a logarithmic scale). With 
the introduction of a vacancy the conductance around 
-E = is dramatically increased to a value of order 
G ~ 0.01(2e^//i). Note that this is the conductance in- 
duced by a single vacancy. In a large system containing 
many such defects one can expect even larger enhance- 
ments of the intertube conductance. If the defects con- 
tribute incoherently (see Sec lVB 3|l . the intertube con- 
ductance will increase proportionately to the number of 
defects. 

We have investigated the dependence on intertube con- 
ductance of the position of the vacancy. (The effect of 
varying the strength of the potential will be investigated 
in the next section.) The inset to Fig.^Jplots the results 
obtained for intertube conductance at i? = by moving 
the vacancy around the circumference of the tube on the 
10 sites of the (10, 10) armchair tube that are crystallo- 
graphically equivalent. The position of the vacancy has 
significant impact on the intertube conductance, with the 
largest enhancement of intertube conductance when the 
vacancy site is closest to the other nanotube. This indi- 
cates that the wave that is backscattered by the vacancy 
in the defective tube has a large amplitude localised in 
the vicinity of the vacancy. 

An experimental investigation of the effects of vacan- 
cies (introduced intentionally by sputtering) on the con- 
ductance of nanotube ropes was carried out by Stahl et 
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FIG. 12: Intershell conductances for DWNTs in the pres- 
ence of a single vacancy (placed on the inner tube), for the 
three types of coupling determined by the angular symmetry 
of jrrni The solid line is a (3,3) to (15,0) system ("zero" 
coupling); the dashed line is a (5, 5) to (18, 0) system ("inter- 
mediate" coupling); and the dotted line is a (7,7) to (21,0) 
( "strong" coupling) . 

aLiii The results were interpreted in terms of the in- 
tertube tunnelling between metallic tubes of matching 
chirality in a rope consisting of a random distribution 
of tubes of random chiralities. Tunnelling between car- 
bon nanotubes of different chiralities was ignored in the 
analysis, based on the existing theorj*^ for clean tubes 
which showed that this tunnelling is heavily suppressed 
as a result of momentum conservation. Our results show 
that, owing to the loss of momentum conservation in 
the vicinity of the vacancy, the intertube transport be- 
tween nearest neighbour tubes of different chirality can 
occur efficiently. Even for a single vacancy, the inter- 
tube conductance between two neighbouring tubes can 
be as large as 0.01(2e^//i), corresponding to a resistance 
of around IMTl. This is well within the r ang e of intertube 
resistances that was measured in Ref. ^l[ Our results 
therefore indicate that owing to the loss of momentum 
conservation at the vacancies, the tunnelling between in- 
commensurate certainly carbon nanotubes cannot be ne- 
glected in these experiments, and is likely to play a sig- 
nificant role. 



2. Multtwall Geometry 

We now consider the effects of a vacancy in a DWNT 
geometry. The vacancy breaks the rotational symmetry 
described in ijllDI thereby removing the distinction be- 
tween the different strengths of coupling of the metallic 
subbands. Therefore, if all other parameters are equal, 
one would expect the differences between "Zero" , "Inter- 
mediate" and "Strong" cases to disappear. Fig. [T^ shows 
the conductance for the same three DWNTs considered 
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in jjlV CI now with a single vacancy in the armchair tube. 
The magnitude of the intertube conductance for all three 
types is now G w 10~^ (2e^//i) around E = 0. Recall- 
ing that the intertube conductances for clean tubes, dis- 
cussed in illV CI differ by many orders of magnitude in 
these three cases, one indeed see that the vacancy does 
break the distinction between these cases giving conduc- 
tances of similar magnitudes. Furthermore, the overall 
size of the conductance is very much increased. Again, 
the presence of even a single vacancy can very signif- 
icantly affect the conductance between the shells of a 
MWNT. 



B. Other Types of Defect 

The nature of the dominant defect scattering in high 
quality nanotubes is not well characterised. Indeed, it is 
likely to vary between systems involving different prepa- 
ration methods and different geometries. The disparity 
between the intershell MWNT conductance experiments 
of vertically suspended MWNTs^ and those where the 
MWNTs are placed on a substratei^ suggests that de- 
fects may arise due to substrate effects. We have studied 
the effects of various types of localised defects, includ- 
ing potential defects and deformations of the outer tube 
arising from substrate roughness. 

We focus on the case of MWNTs, in order to make con- 
nection with the experimental measurements of Bourlon 
et alw^ We have calculated the effects on the conductance 
at = of a single defect in the outer shell of a DWNT, 
consisting of a (13, 13) armchair (outer shell) and a (12, 0) 
zigzag (inner shell) carbon nanotube. According to our 
classification of rotational symmetry, this combination 
is of "intermediate" coupling, which is the typical case. 
The effect of a single defect is determined by choosing a 
sufficiently long section of tube L and sufhciently smooth 
contacts Lc that finite size effects are negligible. In order 
to quantify the effect of the scattering, we report dimen- 
sionless measures of intratube resistance, and intertube 
conductance, defining the quantities 

^. = l-y = 1-^E^m.1;m',2 (22) 

Is = ^ ^ Tfj.,p-,fj.' ,p' (23) 

l,2;p'— 3,4 

Rs is a measure of the intratube backscattering,'^'^ and Is 
is a measure of the intertube conductance. In B 31 we 
shall discuss how these quantities can be used to compare 
with experiment.— 



1. Potential Scatterer 

We consider a potential centred on a site rg of the 
outer carbon nanotube with strength Vq and correlation 



length, ^, such that there is an additional site potential 
on this tube of 

= Foe-l'-'-'-°l'/«' (24) 

This is a model for several types of defect: molecules 
absorbed onto the surface of the MWNT; a residual po- 
tential due to the close proximity to the substrate; ni- 
trogen substitution, which can be modelled by choosing 
Vo w -2.5eV.32 

In Fig. El the dimensionless parameters of Eqns. [23 
and are plotted against potential strength Vq for dif- 
ferent values of the correlation length, ^. In each case, 
the backscattering Rs and intertube coupling Is increase 
at small Vq as Vq consistent with expectations from per- 
turbation theory. In this weak scattering regime, the 
relative sizes of the backscattering and intertube scat- 
tering depend strongly on the correlation length ^. One 
might expect that increasing the range of the potential 
^ at fixed amplitude Vq would cause a larger effect of 
the defect as the total integrated weight of the potential 
increases. This is true for Is and for Rs at large Vq. How- 
ever, for small Vq, the intratube backscattering Rs actu- 
ally falls with increasing ^, even to below Is for ^ — 2.0a: 
in this case, the defect more efficiently causes intertube 
scattering than it does intratube backscattering. This 
strong dependence on the correlation length can be un- 
derstood as a consequence of the suppression of intratube 
backscattering due to the conservation of "pseudo-spin" 
by the (perturbative) scattering from a potential with 
range larger than a lattice constant. ^-^-^ For large ^ the 
backscattering parameter Rs becomes small; there is no 
such suppression of intertube scattering, so Is remains 
large. 

The effects of the potential scatterer are very similar 
for both positive and negative Vq. For ^ = 0, both Rg 
and Is are symmetric in Vq. For ^ 0, an asymmetry in 
Vq arises, but, as shown in Fig. ^1 even for ^ = 2a the 
asymmetry remains very small. 

2. "Dent" Scatterer 

One possible source of scattering in the measurements 
of Bourlon et ali^ is the roughness of the substrate, 
which could lead to deformation of the MWNT.iS, We in- 
vestigate this by considering a "dent" in the outer shell, 
which when placed at a position Vq gives rise to a corre- 
lated displacement of the atoms. We displace the atom 
i (with undeformed position r;) in the radial direction 
towards the tube axis by a distance 

,Sr, =doe"l'''-''«l'/«'. (25) 

These displacements will change the inter- and intra-tube 
hopping parameters and could also lead to additional on- 
site potentials. Provided the curvature of the tube is 
negligible, the correction to the intratube hopping and 
onsite potential is of order (dg/a)^ (by symmetry there 
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FIG. 13: (Color online) The dimensionless measures of 
backscattering, Rs, (solid line) and intertube scattering, Is, 
(dashed line) for a (13,13) to (12,0) DWNT with a single po- 
tential defect, of varying strengths Vb and correlation length 
^ [see Eqn. 1241 . The total length of the system is L = 87a 
with tunnelling introduced smoothly via a contact length of 
Lc = 2a, so boundary effects are minimal and the results are 
characteristic of a single defect in the bulk. 
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FIG. 14: (Color online) The dimensionless measures of 
backscattering, Ka, (solid lines) and intertube scattering. Is, 
(dashed lines) for a (13,13) to (12,0) DWNT with a single 
"dent" in the outer shell as a function of the displacement do 
and correlation length ^ [see Eqn. 1251 . Results are shown for 
the set of correlation lengths ^ = 0.00, 1.76, 3.53, 5.29A, which 
may be assigned to the graphs by noting that both Ks and Is 
increase with increasing ^ for do ~ lA. 



can be no term linear in d^). However, the correction to 
the intertube hopping is of order d^j a. We therefore con- 
sider only the changes to the intertube hopping, which is 
accurate for small deformations do ^ a. 

Fig. 1 141 shows the results of calculations of Rg and I^, 
Eqns. 123 and 1231 for a single dent scatterer of various 
correlation lengths ^ as a function of the strength Vq. 
Compared to the potential scatterers, the "dent" scatter- 
ers give qualitatively different results. Now, the intertube 
conductivity is always higher than the intratube resistiv- 
ity. This is a consequence of the displacement having no 
effect on the intratube tight-binding parameters (valid 
for do <C a): intratube backscattering involves a second 
order process in which the electron tunnels to and from 
the other tube. Thus our results indicate that the dent 
scatterer especially enhances the intertube conductance 
as compared to the intratube resistance. 



3. Experimental Comparison 



In the experimental studies of Bourlon et al. ; non- 
local voltages induced between contacts to the outer 
shell of MWNTs were used to extract the intrashell 
(6 — 25) kfi/im"^, for each of the 



resistivity, pr^l^ 



outer two shells and an intershell conductivity, gf^ter = 
[(3.7 — 20) kfi]" between these two outer shells. 

(Contributions from other shells were shown to be small.) 
As emphasised in Ref. IT^ the intertube conductivity 
cannot be understood within existing theoretical de- 
scriptions of intershell coupling in incommensurate nan- 



B Other Types of Defect 
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otubes. We believe, as suggested in Ref.lT^ that the dis- 
crepancy is due to the effects of disorder which was not 
considered in previous theoretical works. We can use our 
results on impurity scattering to estimate the intrashell 
resistivity and intcrshell conductivity in the experiments. 
Since the nature and density of defects are very uncertain 
in these experimental systems, our approach is to search 
for parameter regimes that are consistent with the exper- 
imental results. This provides constraints on the nature 
of the defects in these systems. 

Comparing the results of our calculations with these 
experimental results^^ requires a number of simplifying 
assumptions. 

While the MWNTs used in the experiments contained 
around 20 shells, since transport was shown to be largely 
confined to the outer two shells we consider just two coax- 
ial shells. We use the results shown above on the (13, 13) 
to (10, 0) system, which is a large radius system showing 
the effects of incommensurability, and with typical "inter- 
mediate" coupling. The overall diameter of this system, 
of about d = 16A, is much smaller than the diameter 
D = 17nm of the outer shell of the MWNTs in the ex- 
periment. We take account of this difference by scaling 
the scattering rates we obtain in the numerics on a sys- 
tem with diameter d by a factor of (d/D)"^ ^ 1/100. S 
[This reduction reflects the reduction in the matrix ele- 
ment of a localised potential as l/D, and thus a reduction 
in scattering rate (using Fermi's golden rule) as l/D"^.] 

In the experimental system the tubes are doped such 
that several modes are occupied, with the estimated num- 
ber of modes A^^^ = 10 - 20M There are then 2N^ Fermi 
points in each shell. Generically, the Fermi momenta of 
adjacent shells will not coincide and the shells will be in- 
commensurate. We take the results of the above calcula- 
tion of the effects of defects on the transport between the 
metallic bands of armchair and zigzag tubes as a measure 
of the scattering probabilities between states at the indi- 
vidual Fermi points of the two outer shells of the MWNTs 
in Ref. We consider the limit of weak backscattering 
and weak intershell scattering from a single defect. In 
our numerical results, each shell has 4 Fermi points, and 
we use the total intertube scattering probability Is to 
obtain the intertube scattering probability between indi- 
vidual Fermi points as /s/4^ (assuming that the 4 Fermi 
points in one tube scatter equally to the 4 Fermi points 
in the other). Then, for N^^ modes, we extend to 2iVin 
Fermi points in total to obtain that the scattering prob- 
ability for a single incoming mode into any of the modes 
on the second tube is 

P/ = (2iV,„)^(-^)' (26) 

which incorporates the reduction of (d/D)^ discussed 
above. Similarly for the reflection probability, extend- 
ing from the case (in the numerics) of two backscattering 
modes to N^i backscattering modes, the mean backscat- 
tering probability for a single incoming mode into any of 



the modes on the second tube is 

Pfi = (^m)^(;^)'. (27) 

We have ignored current conservation in these estimates, 
which would involve the introduction of compensation 
between intertube and backscattering probabilities; this 
assumption is valid provided scattering is weak, and both 
Pi and pr are small. 

The experiments involve the scattering by many de- 
fects, while the above numerical results are for a single 
scatterer. We extend to larger numbers of scatterers by 
assuming that the scatterers combine incoherently, with 
Ug defects per unit length. This will give a good approx- 
imation to the transport properties provided the system 
is not in a strongly localised regime. The experimen- 
tal samples are far from the strongly localised regime 
even down to low temperatures: from the experimental 
measurements the conductance on the lengthscale of the 
contact separation (< Ifim) remains large compared to 
e^h. 

Within these assumptions, and using the analysis pre- 
sented in Appendix IIXI we find that, for rig scatterers 
per unit length, the intertube conductivity gintor and in- 
tratube resistivity pintra (both per unit length) are 



^7inter — 



Pintia — 




Note that, within this weak scattering approximation, 
the resistivity is independent of the number of modes: 
as Nm increases, the increased number of Fermi points 
into which backscattering can occur (increasing the re- 
sistivity) is compensated by the increase in the number 
of modes available for conduction (increasing the con- 
ductance). On the other hand the intertube conductance 
increases as the square of the number of modes, N,^: one 
factor coming from the number of incoming modes, and 
one factor from the increased number of Fermi points into 
which any one mode can be scattered in the other shell. 

We are now in a position to compare our results with 
the experimental resultsii^ It is helpful to note that 
within our theory the ratio of gintor to pintra is indepen- 
dent of the defect density 

/Ointra \ " / JrCs 

In the experimentSfi^ this ratio is of order 5°ntor/ftntra ~ 
(2e^//i) . Agreement of the theory H30|l with this experi- 
mental value requires Rs/Is ^ 0-5N^ ~ (50—200), using 
the estimatei^ A^m ^ 10 - 20. Within the wide range of 
uncertainty in the parameters, the experimental results 
are consistent with our results for potential defects with 
vanishing correlation length ^ = 0, for which we find 
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Rs/Is — 180 over the range of Vq shown in Fig tT^ a). 
or for potentials with a small correlation length. In our 
calculations, the ratio Rs/Is becomes small for potential 
scatterers with large correlation lengths; this reduction 
is likely to be overestimated in a calculation based on 
metallic bands alone. Note that in the case of "dent" 
defects we find that Rs/Is ^ 1 in all cases fFig ll4|l . 
Thus, our analysis shows that the experimental results 
are consistent with the dominant disorder being poten- 
tial scatterers with a small correlation length, but are not 
consistent with the dominant disorder being "dents" in 
the outer tube caused by substrate roughness 

To estimate the strength of the disorder, we focus on 
the case of potential defects with vanishing correlation 
length, ^ = 0. From the results in Fig ll3f al. we find that 
Rs ~ 1.3 X lO""' (Vb/eV) over the range of Vq shown. 
Matching (^i^ to p°^tra ^ lOkil/^m^^ sets an estimate for 
the required number of potential scatterers rig per unit 
length as a function of their strength Vq. Converting to 
the number of scatterers per atom, V3a^rts/(47r_D), we 
find agreement for a root mean square potential energy 
fluctuation on each site of about 0.8eV. This is com- 
parable to the disorder strengths required by other the- 
oretical studies to account for a resistivity of this size. 
[Indeed, Eqn. [SHI closely matches the results of Ref. IsM 
for Nm = 2, relating mean-free-path to resistivity by the 
Drude formula.] 

In summary, our theory shows that the dramatic in- 
crease of intershell coupling arising from scattering from 
disorder is sufficient to account for the experimental ob- 
servations. We find quantitative agreement (within the 
uncertainties of the experimental parameters) for a model 
of disorder consisting of short-range potential scatter- 
ers. We note that our theory shows that the relative 
sizes of the intershell scattering and intrashell backscat- 
tering depend strongly on the number of modes A^m, with 
5inter/pintra oc N"^. It would bc interesting to test this 
prediction in experiment. 



VI. CONCLUSIONS 

We have studied the transport between aligned carbon 
nanotubes (or nanotube shells) in both rope and mul- 
tiwall geometries. We focused on cases of two metallic 
nanotubes of different chiralities, with mismatched Fermi 
momenta and incommensurate periodicities. Our numer- 
ical calculations show clear evidence of the importance of 
momentum conservation, and of a classification we intro- 
duced for the rotational symmetry of a pair of shells in a 
MWNT. For clean tubes, we find that the loss of momen- 
tum conservation at the contacts dominates the intertube 
transport. 

We studied in detail the effects of localised defects 
on the intertube transport in incommensurate nanotube 
ropes and MWNTs. Defects break momentum conserva- 
tion, and the rotational symmetry in multiwall geome- 
tries. We find that the intertube conductance of incom- 



mensurate nanotubes is very dramatically increased by 
even a single vacancy. Our results show that the in- 
troduction of vacancies to nanotube ropes enhances the 
tunnelling between incommensurate tubes sufficiently to 
make this defect-mediated tunnelling an important trans- 
port mechanism. We calculated the effects of localised 
defects of various types on the intra- and inter-shell trans- 
port of MWNTs. For a potential scatterer, we find a 
strong dependence on the range, ^, of the potential: for 
large ^, intratube backscattering is suppressedSiM but in- 
tertube scattering remains large. For a "dent" scatterer, 
involving a localised depression of the outer shell, inter- 
tube conductance is very much enhanced with intratube 
backscattering remaining small. 

Our results establish the extreme sensitivity to disor- 
der of the intertube transport of incommensurate nan- 
otubes with mismatched Fermi momenta. They provide 
the first quantitative description of the effects of this scat- 
tering on the transport between nanotube shells in rope 
and multiwall geometries. A comparison of the results 
of our calculations with experimental measurementsiS 
shows agreement with our theory for a model of disor- 
der consisting of short-range potential scatterers. Poten- 
tial scatterers with large correlation lengths and defor- 
mations of the outer shell are not consistent with the 
experimental results. 

Finally, we would like to comment briefiy on the effects 
of electron-electron interactions. In a one-dimensional 
metal, interactions lead to a smearing of the Fermi sur- 
face, so that there is no longer a discontinuity in the 
occupation numbers at the Fermi surface. Owing to 
the central role played by momentum conservation in our 
considerations, it is important to ask whether this smear- 
ing of the Fermi surface could affect our results. We find 
that including interactions within a Luttinger liquid anal- 
ysis does not change the essential behaviour described 
above for non-interacting electrons!^ the intertube con- 
ductance is still strongly suppressed for incommensurate 
nanotubes with mismatched Fermi wavevectors, and is 
strongly enhanced by impurity scattering. The interac- 
tions will, in addition, cause the overall intertube con- 
ductance to be suppressed as a power law in tempera- 
ture (or bias), as is familiar from tunnelling into a Lut- 
tinger liquid. These interaction effects will be weaker at 
high temperatures and for MWNTs with many occupied 
subbandsi^ 
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VII. GREEN'S FUNCTIONS FOR 
SEMI-INFINITE NANOTUBES 

To determine the Green's functions we make use of the 
form 



E - en,k + iS 



(31) 



where l^"^,*;) are energy eigenstates, /i is the mode index 
and k is the wavevector along the axis of the carbon 
nanotube. Positions are labeUed by i,a, with i labeUing 
the ID unit ceU (along the tube axis), and a the atoms 
within this ID cell. 



A. Armchair Tube 

The metallic bands in an armchair carbon nan- 
otube have energy dispersion relations 



1 



2 cos ( — 



(32) 



where /i = ± labels the two metallic bands, k is wavevec- 
tor along the axis of the nanotube. For a tube with pe- 
riodic boundary conditions, we write the wavefunction 
|^/i,fe) in the form 



(«,a|*M,fe) 



NlD 

2 



(33) 



where Xi is the distance of the ID unit cell along the tube 
axis, iViD is the number of ID unit cells. The functions 
X^,fe are easily obtained from the free-particle eigenstates; 
explicitly for an (n, n) tube they are 



X±,k [a) 



1 

An 



I 1 \ 

±e'^ 

■ka 

— e 2 
V : ) 



(34) 



where the atoms in the ID unit cell are labelled in order 
of their angular positions around the tube. 

For a finite length tube, the wavefunction is a stand- 
ing wave, chosen to satisfy the boundary conditions at 
the two ends of the carbon nanotube, which places re- 
strictions on the allowed values of k. The armchair finite 
tube wavefunctions take the form 



C 



2 e 



X^,fc(a)-e'^e-'=-'X;,. (a) 
(35) 

where C ~ \/8nN is a normalisation factor, which is 
independent of k for a long tube. 

Substituting the wavefunctions and energy dispersion 
relation into Eqn. 1311 we can convert the summation to 
an integral, and using symmetry relations and the fact 
that the resultant integral is analytic at fc = we find 



271- 



d(ka) 



27r 



E 



I 



E - e^,fc -I- i6 



(36) 



The calculation simplifies further by considering only the 
Green's function at the end of the carbon nanotube, i.e. 
i = j = 1. (This is all that is required in our calculations 
of the influence of the semi- infinite leads in The de- 

pendency on the inner product in Eqn. 1361 fdescribing the 
coupling between the different sites labelled by a) can be 
reduced to a 2 x 2 matrix, since there are only two types 
of atom on the "edge" of the armchair nanotube. Using 
contour integration in the variable z — , we find the 
Green's function for an armchair carbon nanotubc"^^ 



\G\d) = 



1 f-z+ + Z- Z+ + Z- 



2nt 



Z+ + Z- -Zj, 



where 



iTE/t 



iTE/t 



(37) 



(38) 



B. Zigzag Green's Functions 

The two metallic bands in an (n, 0) zigzag tube have 
modes, /i, which are distinguished by the subband indices 
= 2n/3,4n/3. They are degenerate with energies 



±2tsin 



(39) 



for the right and left moving modes (±). For a tube 
of finite length, the boundaries induce reflections but no 
mixing between the bands of different q since rotational 
symmetry is preserved. The finite tube wavefunctions 
take the form 



c 



3''"-X^,fe («) 



-ikxi 



-fe(«)] 

(40) 



B Zigzag Green's Functions 
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with the same notation as in the armchair case (but note 
that the ID unit cell and associated wavefunctions x are 
now those appropriate for the zigzag tube) . 

We proceed as in the armchair case by restricting our 
attention to the relevant atoms on one end of the carbon 



nanotube. Converting the summation to an integral and 
using residue calculus we find the total metallic Green's 
function on the end of a semif-infinite zigzag carbon nan- 
otube to be 



(a\G\ 



nt 



E 
2t 



2i, 



(41) 



r 



where the terms of two periodicities in j/q — arise from 
the An/ 3 and 2n/3 bands. 



VIII. TRANSMISSION PROBABILITIES 



nanotubes. Following HllU the input wave is labelled by 
mode fj, and lead p, and the output by mode fi' and lead 
p'. One must distinguish between the cases where the 
input and output leads are of the differing chiralities 



Here we state the transmission probabilities across a 
scattering region between armchair and zigzag carbon 



J 



p armchair, p' armchair Ti^.p-fj.',p' 



p zigzag, p' zigzag Tf,^p.^>^p> 



p armchair, p' zigzag (or vice versa) 7)t,p;/^',p' = 



2h 

at V V3 



(42) 
(43) 

(44) 



We use the terminology and definitions for the carbon 
nanotube wavefunctions described in Appendix IVIII and 
Uac and Vzz are the velocities for armchair and zigzag 
tubes. The inner product of the Green's function with 
the ID wavefunctions XtJ. X/j' are to be taken on any ID 
unit cell in the leads p and p' (wavevector labels are 
suppressed, but should be understood to be the Fermi 
wavevectors of the leads). The Kronecker delta func- 
tions represent the modification required for reflection 
back into the same lead and mode. 



IX. INCOHERENT SCATTERING 
SUMMATION 

We illustrate the method used to combine probabil- 
ity s-matrices, which enables us to calculate conductance 
through many scattering centres incoherently, as used in 
i|VB3l Consider the setup in Fig. E| For two identical 
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FIG. 15; Combining two identical scattering sections 



scattering sections we have 

(i:) - (TS)(ij:) («) 
(^:)^(?S)(tr) <«) 

where Lout is a vector representing the probability flux 
of all the outgoing modes to the left of the scatterer, and 
Lin is a vector representing the probability flux of all the 
ingoing modes to the left, similarly for M and R. R and 
T are matrices representing the reflection and transmis- 
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sion of different leads. We assume that the scattering 
sections are symmetric in space. We can eliminate the 
M terms to obtain the combined structure 



Lout 
Rout 



R' T' \ / Li„ 
T' R' j I 



where 



R' 



T' = T [I 
R + TR [I 



R^] 
R^l 



(47) 

(48) 
(49) 



For the setup with two carbon nanotubes, it is useful to 
separate tubes A and B: 
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(50) 



where I is the intertube transmission probability, which 
has been set equal between all modes to simplify the 



problem. Taking one mode per tube, adding TVs iden- 
tical scatterers, and expanding in a power series in terms 
of /, the intertube transmission probability, the leading 
behaviour of the relevant probabilities is 



[Ns ~ 1)T 



lim 



T 



oo Ns{l ~ T) 



(51) 
(52) 



where T^^ is the total intratube transmission probability 
through TVs scatterers if T is the probability of transmis- 
sion through one scatterer, and In^ is the total intertube 
transmission probability through Ng scatterers if / is the 
probability of transmission through one scatterer. We 
make use of these results in 3V B 31 in the limit of weak 
scattering (1 - T), / < 1. 
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